超声速膨胀角入射激波/湍流边界层干扰直接数值模拟

为了揭示膨胀效应对激波/湍流边界层干扰区内复杂流动现象的影响规律,采用直接数值模拟方法对来流马赫数2.9、30°激波角的入射激波与10°膨胀角湍流边界层相互作用问题进行了数值研究。系统地探讨了激波入射点分别位于膨胀角上游、膨胀角角点和膨胀角下游3种工况下膨胀角干扰区内若干基本流动现象,如分离泡、物面压力脉动及激波非定常运动、湍流边界层统计特性和相干结构动力学过程等。结果表明,激波入射点流向位置改变对分离区流向和法向尺度的影响显著,尤其是当激波入射点位于角点及其下游区域。研究发现,膨胀角干扰区内物面压力脉动强度急剧减小,分离区内压力波向下游传播速度将降低而在膨胀区内将升高,膨胀效应极大地抑制了分离激波的低频振荡运动。相较于入射激波与平板湍流边界层干扰,入射激波流向位置改变对膨胀角再附区速度剖面对数区及尾迹区影响显著,将导致其内层结构参数升高而外层降低,近壁区内将呈现远离一组元湍流状态的趋势。此外,流向速度脉动场本征正交分解分析指出,主模态空间结构集中在分离激波及剪切层根部附近而高阶模态以边界层内小尺度正负交替脉动结构为主。低阶重构流场结果表明,前者对应为分离泡低频膨胀/收缩过程而后者表征为分离泡高频脉动。     

基于RANS/LES混合方法的分离流动模拟

飞行器在大迎角、快速俯仰机动时,流场中含有大尺度、非定常的涡结构,传统雷诺平均Navier-Stokes （RANS）模型不能准确模拟流场结构,根据国际上相关研究的发展趋势,需要采用混合RANS/大涡模拟（LES）模型来对复杂分离流动进行准确模拟。本文对基于分区混合与湍流尺度混合的双重RANS/LES混合计算模型进行发展与应用。通过典型简化模型的静、动态湍流大分离流动,测试和验证所采用的脱体涡模拟（DES）类方法,重点研究改进的延迟DES （IDDES）模型在动态问题应用中的正确性和有效性,并对所采用的数值模拟方法和相应的计算软件的可靠性、鲁棒性以及精度进行了考核验证。典型算例包括超声速圆柱底部流动、跨声速方腔流动、NACA0015机翼深失速分离涡模拟等。计算表明：发展的IDDES类混合计算模型可有效解决对数层不匹配的问题;对于定态非定常分离流动,DES、DDES、IDDES等模型计算结果差别不大,随着流动的非定常特性增强,IDDES模型的优势逐渐显现;对于动态非定常分离流动,则需要采用IDDES类模型。     

三维自适应终端滑模协同制导律

针对多枚导弹协同作战的问题,且多枚导弹之间保持有向拓扑通信的条件下,基于终端滑模法设计了视线方向及视线法向的双层协同制导律。其中,视线方向的制导指令能够保证导弹同时完成拦截任务;视线法向上的三维制导律能够保证每枚导弹以期望的视线角攻击目标,从而发挥各枚导弹的最大杀伤力,并且视线角的约束相当于规划了末制导段导弹的弹道问题,在一定程度上避免攻击目标前导弹间发生碰撞。同时,针对所设计的滑模制导律设计了新的自适应律,从而加快了滑模面的收敛速度并且削弱了由符号函数引起的系统抖振现象。基于李雅普诺夫稳定性理论,证明了所设计制导律的正确性,并在最后给出了数学仿真实验,验证了所设计制导律的有效性及优越性。     

膨胀效应对激波/湍流边界层干扰的影响

采用直接数值模拟方法对来流马赫数2.9,30°激波角的入射激波与膨胀角湍流边界层干扰问题进行了数值研究。入射激波在壁面上的名义入射点固定在膨胀角角点,膨胀角角度分别取为0°、2°、5°和10°。通过改变膨胀角角度,考察了膨胀效应对干扰区内复杂流动现象的影响规律和作用机制,如分离泡、物面压力脉动特性、膨胀区湍流边界层和物面剪切应力脉动场等。研究发现,膨胀角角度的增大使得分离区流向长度和法向高度急剧降低,尤其是在强膨胀效应下分离泡形态呈现整体往下游偏移的双峰结构。物面压力脉动功率谱结果表明,膨胀角为2°和5°时,分离激波的非定常运动仍表征为大尺度低频振荡,而膨胀角为10°,强膨胀效应极大地抑制了分离激波的低频振荡,加速了下游再附边界层物面压力脉动的恢复过程。膨胀区湍流边界层雷诺剪切应力各象限事件贡献和出现概率呈现逐步恢复到上游湍流边界层的趋势,Görtler-like流向涡结构展向和法向尺度变化剧烈,同时在近壁区将诱导生成大量小尺度流向涡。此外,物面剪切应力脉动场的本征正交分解分析指出,膨胀效应的影响体现在低阶模态能量的急剧降低从而使得高阶模态的总体贡献相对升高。     

变循环发动机过渡态性能直接模拟方法

为了实现多变量可调的变循环发动机（VCE）过渡态性能模拟和控制规律设计,在稳态逆算法和过渡态性能隐式格式计算方法的基础上开发了变循环发动机过渡态直接模拟方法。在稳态逆算法模型中加入容积效应和转子动力学方程,实现了同时给定加速率、涡轮前温度和压缩部件工作点的条件下直接模拟计算过渡态过程中变几何参数和燃油流量的调节规律,验证了方法的精度和可行性,并将该方法用于变循环发动机转模态性能模拟和控制规律设计。计算结果表明,过渡态直接模拟方法的误差在0.58%以内,超声速巡航状态下由单外涵模态转换到双外涵模态的时间约为1 s,海平面静止状态下双外涵转单外涵的时间约为3 s,且推力、转速、涡轮前温度、喘振裕度等参数过渡平稳。该方法可简化转模态控制规律设计流程,并提高设计精度。     

0.3 m低温风洞液氮供给系统研制

基于系统级的一维热流体模拟分析,建立了适用于研究分析0.3 m低温风洞液氮供给系统的数学模型,并开展了系统漏热、两相流及缓冲罐中液氮容积等流体动力学分析;在系统现有控制策略及试验数据的基础上,基于该数学模型开展了系统压力动态响应分析,获得了在阀门动态调节过程中管网压力的瞬态响应,计算结果与试验值的总体误差控制在10%以内。喷射压力一致化改造避免了阀间干扰,添加的回流管道消除了供给末端的两相流现象,使喷射压力控制精度达到1.1%,调节时间减少到23 s,实现了风洞总温快速安全调节和精确控制。     

非正交弧线齿面齿轮齿面设计及根切研究

非正交弧线齿面齿轮是以端面渐开线弧线齿圆柱齿轮为假想刀具包络展成的面齿轮。根据面齿轮啮合原理推导了非正交弧线齿面齿轮齿面方程,结合Matlab数学模型和根切理论研究了该种面齿轮的根切现象,计算出面齿轮不产生根切的最大内径,并通过Catia仿真滚齿来验证计算的准确性。该研究旨在提出一种新型非正交面齿轮,通过理论推导及仿真模拟来研究该种面齿轮的根切现象,在此基础上探究非正交弧线齿面齿轮的内径影响因素。研究结果表明:在同等参数下,当位置参数增大时,刀具渐开线截面高度和内径逐渐减小,当刀具圆弧半径增大时,内径逐渐减小,刀具渐开线截面高度逐渐增大。     

导流槽构型对火箭起飞时流场和声场的影响

为了研究四喷管运载火箭起飞时火箭周围噪声环境问题,建立燃气/空气双组分的可压缩流动模型,采用2阶Roe格式、SAS(scale-adaptive simulation)湍流模型和声学类比积分法Ffowcs-Williams Hawkings(FW-H)求解三维Navier-Stokes方程。以单机火箭的噪声问题为对象开展数值模拟,并将噪声数值计算结果与试验数据对比,误差在3dB以内(相对误差小于1.6%),验证预测噪声方法的有效性,进而研究四喷管运载火箭起飞阶段采用不同导流槽构型对箭体舱段区域噪声环境的影响,结果表明:在相同噪声接收点处,单侧导流槽对应的总声压级比双侧导流槽大,两者的总声压级(OASPL)之差最大为10.7dB。另外,当采用单侧导流槽时,沿周向接收点的噪声总声压相对单侧导流槽中心截面呈对称分布,而且沿导流出口方向逐渐增大。所建立的噪声数值方法为大推力捆绑运载火箭舱段的噪声环境预测及其控制提供一定参考。     

发动机唇口电热防冰系统性能仿真

采用三维耦合传热方法对发动机唇口电热防冰系统进行了多状态下的性能仿真,获得了唇口表面平衡温度及动态温升响应时间,以此评估防冰系统性能,并通过与冰风洞试验结果对比,对仿真过程进行了修正。研究表明:唇口表面加热区平衡温度均高于273.15 K,温升响应时间小于60 s,满足防冰需要,较为严酷的防冰状态点主要为飞行速度较大、环境温度较低的状态点;修正唇口材料导热系数后,仿真与试验结果吻合良好。研究结果能有效指导工程研制,对提高防除冰仿真分析结果的准确性提供了很好的借鉴。     

退化量有缺失的无失效小样本轴承可靠性评估

针对有试验性能退化量参数记录、且(或)有部分数据缺失情况的小样本无失效轴承试验问题,通过由Taylor和Thompson提出的数据模拟法实现补全样本退化量,结合Bootstrap自助法扩大样本量,再根据基于性能退化轨迹的补充信息方法来进行其可靠性评估。选取7组受试轴承的振动退化量,对比在完整数据和带有缺失数据情况下的分析结果,发现可合理利用原舍弃不合规试验的部分有效信息,使之增加可靠性评估的样本数,从而得到较不用这些数据更为准确的结果,且所得结果较用完整数据结果绝对值相差在01以内。对比由极大似然估计法和加权E-Bayes法分析试验寿命数据的结果,发现该方法所得评估结果更优,与试验实际相比误差在10%以内,对于提高评估精度及降低试验成本有积极的实际意义。